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Charged-particle pseudorapidity densities are presented for the d+Au reaction at ^/snn = 
200 GeV with -4.2 < rj < 4.2. The results, from the BRAHMS experiment at RHIC, are shown for 
minimum-bias events and 0-30%, 30-60%, and 60-80% centraiity classes. Models incorporating both 
soft physics and hard, perturbative QCD-based scattering physics agree well with the experimental 
results. The data do not support predictions based on strong-coupling, semi-classical QCD. In the 
deuteron-fragmentation region the central 200 GeV data show behavior similar to full-overlap d+Au 
results at y^ijvjv = 19.4 GeV. 
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The saturation of initial parton densities in relativis- 
tic heavy-ion collisions, a manifestation of high-density 
QCD, is expected to significantly influence the pseudora- 
pidity and centraiity dependence of the emitted charged- 
particle densities from these reactions 0, S ^IS- 
Earlier charged-particle pseudorapiditv density distribu- 
tions for Au-HAu colhsions 7, s', H E from the 
Relativistic Heavy-Ion Collider (RHIC) have been used 
to constrain model predictions for ultrarelativistic heavy- 
ion collisions. They have been inconclusive, however, 
as to whether parton saturation in the initial state con- 
tributes significantly to the reaction dynamics, with both 
saturation-model 0, 0| and calculations that instead fo- 
cus on the energy-loss mechanisms for the multiple mini- 
jets created in the collisions [l^, ll^ llJ, successfully 
describing the data. A similar model ambiguity found 
in explaining the observed suppression of high-pt par- 
ticles in Au-|-Au collisions was recently resolved with 
mid-rapidity d-|-Au data sho wing the suppression is not 
an initial-state effect 0, 0, 0, ^| . It has been sug- 
gested sj that global particle yields in d-f Au collisions 
might result in a definitive signature of parton saturation. 

We report on a measurement of the charged-particle 
pseudorapidity densities for the d-f Au reaction at 
y/SNN = 200 GeV with pseudorapidity 77 coverage of 
—4.2 < 77 < 4.2. The pseudorapidity densities are 
reported for minimum-bias events and 0-30%, 30-60%, 
and 60-80% centraiity classes. The results allow for 



a detailed comparison to model predictions of particle 
production at RHIC energies. The most central data 
(0-30%), where both deuteron nucleons are expected 
to participate in the reaction, are compared with full- 
overlap d-|-Au data obtained by the NA35 collaboration 
at y/s^N^ 19.4 GeV 

The present analysis employs several of the BRAHMS 
global detector subsystems: The Si Multiplicity Ar- 
ray (SiMA) and the scintillator Tile Multiplicity Array 
(TMA) are used for centraiity determination and 
to measure the pseudorapidity densities close to mid- 
rapidity. The Beam-Beam Counter (BBC) arrays are 
used to reconstruct the collision vertex and to determine 
the pseudorapidity densities at larger pseudorapidities. 
The "Inelasticity Counters" (INEL), developed for the 
pp2pp experiment ^23| , are used for a close-to- minimum- 
bias experiment trigger and to provide vertex position in- 
formation in cases where the beam-beam counter arrays 
arc not able to establish this information. Full details of 
the BRAHMS apparatus can be found in ref. ^23] . 

The layout of the SiMA and TMA detectors for the 
d-|-Au experiment is similar to that presented for ear- 
lier measurements of Au+Au multiplicities at ,ysNN — 
130 GeV and 200 GeV, and details of the analysis pro- 
cedures can be found in refs. The SiMA was 
configured with 25, 4 cm x 6 cm Si wafers in an hexago- 
nal arrangement around the beam pipe, with each wafer 
functionally divided into 7 discrete segments along the 
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beam line and located 5.3 cm from the beam axis. Four 
sides of the hexagonal array were populated with six de- 
tectors, each, with the remaining two sides left largely un- 
populated except for a single wafer mounted outside the 
acceptance of either of the BRAHMS spectrometers (2^ . 
The TMA was populated with 38, 12 cm x 12 cm plastic 
scintillator tiles with fibre-optic readout located 13.7 cm 
from the beam axis. The hexagonal TMA array had four 
sides fully populated with eight detectors, each, with two 
and four detectors mounted on the other two sides, re- 
spectively. With this arrangement, the SiMA and TMA 
can each cover the pseudorapidity range —2.2 < ij < 2.2 
for collisions at array center. In the analysis, a range 
of collision vertex locations z about the nominal array 
center is used, with —15 cm < z < 15 cm. Particle 
multiplicities were deduced for an individual SiMA and 
TMA element by using GEANT simulations to con- 
vert the observed energy-loss signal to the number of pri- 
mary particles hitting that element. The HIJING event 
simulator |12| was used to obtain the initial distribution 
of particle types and momenta. 

Two Beam-Beam Counter arrays, positioned around 
the beam pipe on either side of the nominal interaction 
point at a distance of 2.20 m, are used to extend coverage 
out to 77 = ±4.2. Each array consists of separate sets of 
of small (19 mm diameter) and large (51 mm diameter) 
Cherenkov UV-transmitting plastic radiators coupled to 
photomultiplier tubes. Leading particles timing achieves 
a vertex position resolution of « 2 cm. Charged-particle 
multiplicities are deduced from the number of particles 
hitting each detector, as found by dividing the measured 
ADC signal by that corresponding to a single incident 
particle. 

Three pairs of INEL Counters were used to develop a 
near-to-minimum bias trigger by detecting charged par- 
ticles in the pseudorapidity range 3.2 < \ri\ < 5.3. The 
basic INEL counter consists of a plastic scintillator ring 
that is segmented into four pieces and arranged about 
the beam pipe. The counter locations with respect to 
the nominal vertex and (inner radius; outer radius) were 
±155 cm (4.13 cm; 12.7 cm), ±416 cm (6.67 cm; 12.7 cm) 
and ±660 cm(6.67 cm; 12.7 cm). Relative time-of-flight 
measurements with the INEL arrays lead to an interac- 
tion vertex determination with a resolution of ~ 5 cm. 
The INEL counters also provide a minimum-bias trigger 
for the experiment. Based on GEANT simulations, they 
are sensitive to 91 ±3% of the total inelastic cross section. 

Reaction centrality is determined using a geometry- 
weighted average of SiMA and TMA multiplicities. Both 
the SiMA and TMA multiplicities are corrected for the 
distance of the actual interaction vertex from the nominal 
vertex at array center. GEANT simulations were used to 
correct for the possibility that neither the SiMA nor the 
TMA detectors will be hit by a particle for the most 
peripheral events. 

Figure n shows the normalized SiMA and TMA aver- 
aged multiplicities. The observed falloff is unlike that 
of the corresponding Au+ Au spectrum 9] , where there is 




Normalized Si + Tile lUlultiplicity 

FIG. 1: SiMA and TMA averaged multiplicity distribution 
normalized to the 1% centrality level. Lines show efficiency 
corrected limits for indicated centralities. The insert shows 
the correlation between the SiMA and TMA multiplicities. 



an extended "flat" region and a well-defined high mul- 
tiplicity knee. The d+Au spectrum instead reflects a 
smaller number of participants, making the measure- 
ments more sensitive to the underlying nucleon-nucleon 
collision multiplicity distribution. Also, the small parti- 
cle multiplicities for d±Au collisions result in a relatively 
large range in the fraction of particles detected for a given 
total number of particles emitted. The broad correlation 
band found in comparing SiMA and TMA multiplicities 
(Fig. n insert) illustrates the statistical scatter of the 
semi- independent multiplicity measurements. Reaction 
centralities are found by integrating the yield under the 
multiplicity curve. Limits for the 30%, 60%, and 80% 
centrality cuts are indicated by the vertical lines. 

Figure |21 shows the resulting charged-particle 
pseudorapidity-density plots for minimum-bias events 
and 0-30%, 30-60%, and 60-80% centrality classes. The 
SiMA and TMA results have been averaged. Overall 
statistical uncertainties are indicated or are smaller than 
the data points. Systematic uncertainties, denoted by 
the horizontal brackets and estimated as 8% for the 
averaged SiMA and TMA results and 12% for the BBC 
values, are determined by exploring the variation of the 
deduced pseudorapidity densities to reasonable changes 
in the energy calibrations and background subtrac- 
tion. Our minimum-bias data agree withiii systematic 
uncertainties with recently reported results |25| . 

Three model calculations are compared to the data. 
The solid curves show the predictions of HIJING [T^ . 
a Monte Carlo model that includes both soft and hard, 
perturbative QCD-based scattering effects. The dashed 
curves show the predictions of the AMPT model 0, 
0, 13 which includes both initial partonic and final 
hadronic interactions. For comparison with experiment, 
both the HIJING and AMPT model results have been fil- 
tered through a GEANT 24] simulation of the BRAHMS 
experimental response. Centrality is based on the frac- 
tion of events with the highest particle multiplicity within 
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FIG. 2: a) and b) Charged-particle pseudorapidity densities 
for indicated centrality ranges, c) Multiplicity ratios R"""^" 
(squares) and (triangles), as discussed in the text. 

Statistical uncertainties are indicated by vertical lines or are 
smaller than the symbols. Detached horizontal brackets indi- 
cate the total (statistical and systematic) uncertainties. The 
solid, dashed and dotted curves in a) and b) are the results 
of the HIJING, AMPT and Saturation models, respectively. 
The curves in c) show the HIJING resuhs for R°"^° (solid) 
and (dashed), with the arrows indicating the values 

expected for Au- and d-participant, only, scaling. In all pan- 
els, the connected open circles (asterisks) correspond to unre- 
stricted HIJING calculations with centrality classes based on 
multiplicity (impact paramter), as discussed in the text. 



the pseudorapidity range of the SiMA and TMA ar- 
rays. Both models reproduce the experimental results 
at midrapidity and at positive rapidities approaching the 
deuteron fragmentation region. At negative rapidity (Au 
fragmentation region) the two models start to diverge 
and here HIJING appears to be in slightly better agree- 
ment vi^ith our results. 

The dotted curves in Fig. |21 show the expectations of 
the Saturation Model 5] which accounts for the high- 
density QCD effects that are expected to limit the num- 
ber of partons in the entrance channel. In this case the 
centrality dependence was based on the published curves 
of charged-particle pseudorapidity densities for different 
centrality ranges given in ref. The model appears 
to be unsuccessful in reproducing either the centrality or 
pseudorapidity dependence of the present results. 
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particle densities are shown in Fig. 12; along with the 
corresponding HIJING ratios (curves). Here we take 
(Npart) = 13.6 ±0.3, 8.5±0.3, and 4.7±0.3 for the 0-30%, 
30-60%, and 60-80% centrality ranges, respectively. The 
systematic uncertainties for the ratios include the partic- 
ipant scaling uncertainty and a 5% uncertainty for the 
experimental pseudorapidity density ratios. The partici- 
pant ratios appropriate for Au(left arrows)- and d(right 
arrows)- participant-only scaling are shown in Fig 2c for 
j^o-30 (goij^j) R30-60 (dashed), respectively. The HI- 
JING model reproduces well the experimental ratios, as 
shown in Fig.|3;. In this regard, it can be noted that the 
midrapidity pseudorapidity densities obtained in a stand- 
alone HIJING calculation scale roughly as the number of 
Au participants. 

We use HIJING/GEANT to explore the potential bias 
introduced by the limited acceptance of the MA (SiMA 
and TMA) on the deduced pseudorapidity distributions 
and R values. The connected open circles in Fig.|2show 
the results for the 0-30% (panel a) and 60-80% (panel 
b) centrality classes of an unrestricted HIJING calcula- 
tion where the centrality is based on all charged particles 
emitted in the reaction, and not just those that satisfy 
the experimental acceptance. The greatest effect on the 
pseudorapidity distributions is found for the 60-80% cen- 
trality cut and amounts to as much as an 18% enhance- 
ment in the measured dN/dyy values in the Au fragmen- 
tation region. The corresponding R^"'^" curve no longer 
shows evidence of the mid-rapidity maximum observed 
for the experimental results. 

In lighter systems, the events selected in a given range 
of multiplicity-based centrality are not all the same as 
would be selected if it were possible to base centrality on 
the impact parameter. This is shown in Fig. |21 where the 
connected asterisks indicate the HIJING model dNch/dij 
distributions for 0-30% (panel a) and 60-80% (panel b) 
centrality classes based on impact parameter, and the 
corresponding R°~'^*' ratio (panel c). Here the R°~'^° 
curve shows a steady rise from the d- to Au- fragmenta- 
tion sides, illustrating that the centrality selection does 
affect the deduced l/Npart x dNch/drj values. 

The d-l-Au system was previously studied by the NA35 
experiment at ^/sWn = 19.4 GeV where data were 
obtained for both negative hadrons /i~and for the net 
baryons as measured by the difference of proton and 
antiproton yields (p — p). Pseudorapidity distributions 
were found corresponding to the most central 43% of 
the total inelastic cross section, with both deuteron nu- 
cleons acting as participants. For comparison with the 
present results, the total charged-particle densities are 
determined for the lower energy data taking dNch/drj — 
2 X dN {h~)/dT] + dN {p — p)/dr]. The pseudorapidity den- 
sities are deduced from quoted rapidity distributions by 
first shifting to the center-of-mass system and then as- 
suming the 7r~ mass for h~ and the observed mean-pf 
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FIG. 3: Comparison of central y^sjvjv = 200 GeV re- 
sults(solid squares) with NA35 data (open squares) at 
y^sjvjv ~ 19.4 GeV in a) the nucleon-nucleon center-of-mass 
system, b) the Au rest frame, c) the deuteron rest frame. 
The sohd curve is based on data for Au+Au 0-30% central 
events at ^s]^ = 200 GeV Npart scaling is applied as 
indicated in panel b). 



values for the h~ and p — p distributions. At the higher 
energy of the current measurement, HIJING model sim- 
ulations indicate that the criteria that both deuteron nu- 
cleons act as participants is well satisfied for the 0-30% 
centrality range. In this case it is interesting to see if 
the limiting fragmentation behavior previously observed 
in Au-hAu yields at ^/sWn = 130 and 200 GeV [lIToll is 
present in the d-|-Au system. 

Figure |3t compares the d-|-Au pseudorapidity distri- 
butions in the nucleon-nucleon center-of-mass system, 
where the fixed-target NA35 results have been shifted 
by the center-of-mass rapidity. A factor of 2.2 increase 
is seen in the charged-particle density at mid-rapidity 
for the higher energy data. Using HIJING to determine 
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the different methods of event selection for the two exper- 
iments might reduce this value, the data do not support 



simple participant scaling with energy. Fig. compares 
the data at the two energies in the frame of the Au frag- 
ment, with the NA35 results scaled up by the ratio of 
Au participants at the two energies. The two distribu- 
tions have similar values approaching the Au rapidity, 
although it should be noted that counting the number 
of participants in the heavier fragmentation region of 
a very mass- asymmetric reaction is difficult because of 
multiple scatterings in the target spectator matter [20l |. 
The solid curve shows the results for the 0-30% central 
Au-|-Au distribution at ^/snn = 200 GeV ^3|> scaled 
by the ratio of d-|-Au gold participants to the number of 
Au-|-Au participant pairs. The current measurements do 
not extend close enough to the beam rapidities to assess 
limiting fragmentation scaling on the Au fragmentation 
side. Fig. I^l^ compares the two d-|-Au distributions in 
the deuteron frame. With the given centrality selections, 
Nparti<^) ~ 2 at both energies and so no participant scal- 
ing is done for the comparison. The two distributions 
are found to overlap from roughly one to two units of 
pseudorapidity below beam rapidity, suggesting a limit- 
ing fragmentation behavior over at least this range. 

Pseudorapidity densities of charged particles for the 
d-l-Au reaction at ^/snn — 200 GeV are presented for 
different centrality ranges. The ratio of particle densi- 
ties for central and peripheral events is found to agree 
well with participant scaling in terms of the respective 
fragments away from mid-rapidity. Overall, model calcu- 
lations based on both soft physics and perturbative QCD 
(HIJING, AMPT) lead to excellent agreement with the 
experimental results. Calculations based on the satura- 
tion picture using scale parameters set by previous exper- 
imental data fail to reproduce the measurements and lead 
to a pseudorapidity dependence very different from that 
observed with the current data. Comparison with lower 
energy d-|-Au data suggests a limiting fragmentation-like 
behavior near the rapidity of the deuteron fragment. 
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